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Abstract 



We study the interplay between the local geometric properties and the non-blowup 
of the 3D incompressible Euler equations. We consider the interaction of two perturbed 
antiparallel vortex tubes using Kerr's initial condition [Phys. Fluids 5 (1993), 1725]. 
l/^ I We use a pseudo-spectral method with resolution up to 1536 x 1024 x 3072 to resolve 

■ the nearly singular behavior of the Euler equations. Our numerical results demonstrate 

that the maximum vorticity does not grow faster than double exponential in time, up to 
\^ I t = 19, beyond the singularity time t = 18.7 predicted by Kerr's computations [TH[T7|. 

' The velocity, the enstrophy and enstrophy production rate remain bounded throughout 

the computations. As the flow evolves, the vortex tubes are flattened severely and 
turned into thin vortex sheets, which roll up subsequently. The vortex lines near the 
I region of the maximum vorticity are relatively straight. This local geometric regularity 

' of vortex lines seems to be responsible for the dynamic depletion of vortex stretching. 



1 Introduction 

One of the most challenging questions in fluid dynamics is whether the incompressible 3D 
Euler equations can develop a finite time singularity from smooth and bounded initial data. 
From a theoretical point of view, the main difficulty is due to the presence of the vortex 
stretching term in the vorticity equation, which is formally quadratic in vorticity. If such 
quadratic nonlinearity persists in time long enough, we would expect a finite time singularity 
of the form 0{T — t)~^ in vorticity. Such blow-up rate is consistent with the well-known 
result of Beale-Kato-Majda |2| (see also fl]). There have been many computational efforts 
in searching for finite time singularities of the 3D Euler and Navier-Stokes equations, see e.g. 
[il221IIHllinil23IIlllSlEliniliaiIIll211- One example that has been studied extensively 
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is the interaction of two perturbed antiparallel vortex tubes. This example is interesting 
because of the vortex reconnection which has been observed for the corresponding Navier- 
Stokes equations. It is natural to ask whether the 3D Euler equations would develop a finite 
time singularity in the limit of vanishing viscosity. 

In 113] , Kerr presented numerical evidence which suggests a finite time singularity of the 
3D Euler equations for two perturbed antiparallel vortex tubes. In Kerr's computations, he 
used a pseudo-spectral discretization in the x and y directions, and a Chebyshev discretiza- 
tion in the z direction with resolution of order 512 x 256 x 192. His computations showed that 
the growth of the peak vorticity, the peak axial strain, and the enstrophy production obey 
(T — t)^^ with T = 18.9. Self-similar development and equal rates of collapse in all three 
directions were shown (see the abstract of JH)- While velocity blowup was not documented 
in [131, Kerr showed in his subsequent papers [121 UHl UII that velocity field blows up like 
0{T — with T being revised to T = 18.7. Kerr's computations have generated a lot 

of interests and his proposed initial conditions have been considered as "the most attractive 
candidates for potential singular behavior" of the 3D Euler equations (see page 187 of [T?l]). 

Vortex reconnection of two perturbed antiparallel vortex tubes has been studied exten- 
sively in the literature. Substantial core deformation has been observed [221 HI 131 UHl EUl ESj • 
Most studies indicated only exponential growth in the maximum vorticity. However, the 
work of Kerr and Hussain in fTHl suggested a finite time blow-up in the infinite Reynolds 
number limit, which motivated Kerr's Euler computations mentioned above. 

There has been some interesting development in the theoretical understanding of the 3D 
incompressible Euler equations. It has been shown that the local geometric regularity of 
vortex lines can play an important role in depleting nonlinear vortex stretching f7\ IR| l^t ITIIj. 
In particular, the recent results obtained by Deng, Hou, and Yu [HI HO] show that geometric 
regularity of vortex lines, even in an extremely localized region containing the maximum 
vorticity, can lead to depletion of nonlinear vortex stretching, thus avoiding finite time 
singularity formation of the 3D Euler equations. To obtain these results, Deng-Hou-Yu 
[HI El explored the connection between the stretching of local vortex lines and the growth 
of vorticity. In particular, they showed that if the vortex lines near the region of maximum 
vorticity satisfy some local geometric regularity conditions and the maximum velocity field 
is integrable in time, then no finite time blow-up is possible. See Section 4.2 for the detailed 
description of these results. Kerr's computations fall in the critical case of the non-blowup 
theory in [HI iHIl. To get a definite answer in this critical case, we need to check whether 
certain scaling constants, which describe the local geometric properties of the vortex lines, 
satisfy an algebraic inequality. However, such scaling constants are not available in pHj . 
This is our original motivation to repeat Kerr's computations. 

It is worth mentioning that the predicted singularity time in Kerr's computations is 
T = 18.7, while his computations from t = 17 to t = 18, as mentioned in jj^, seem to be 
under-resolved and were "not part of the primary evidence for a singularity". Clearly, the 
computations for t < 17, which Kerr used as the primary evidence for a singularity, is still far 
from the predicted singularity time, T = 18.7. In order to justify the predicted asymptotic 
behavior of vorticity and velocity blowup, one needs to perform well-resolved computations 
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much closer to the predicted singularity time. Such well-resolved computations can also 
provide more accurate geometric properties of vortex lines, which can be used to check 
whether the non-blowup conditions in |H1 CHI are satisfied. 

In this paper, we perform well-resolved computations of the 3D incompressible Euler 
equations using the same initial condition as the one used by Kerr in [14j. A pseudo-spectral 
method with a very high order Fourier smoothing is used to discretise the 3D incompressible 
Euler equations in all three directions. The time integration is performed using the classical 
fourth order Runge-Kutta method with adaptive time stepping to satisfy the CFL stability 
condition. We perform a careful numerical study to show that the pseudo-spectral method 
we use provides more accurate approximations to the 3D Euler equations than the pseudo- 
spectral method that uses the standard 2/3 dealiasing rule. We use up to 1536 x 1024 x 3072 
space resolution in order to resolve the nearly singular behavior of the 3D Euler equations. 

Our numerical results demonstrate that the maximum vorticity does not grow faster than 
double exponential in time, up to t = 19, beyond the singularity time t = 18.7 predicted 
by Kerr's computations ^[Ej. There are three distinguished stages of vorticity growth in 
time. In the early stage for < t < 12, the maximum vorticity grows only exponentially 
in time. During the intermediate stage for 12 < t < 17, the two vortex tubes experience 
tremendous core deformation and become severely fiattened. Each vortex tube effectively 
turns into a vortex sheet with rapidly decreasing thickness. During this stage, the maximum 
vorticity grows slightly slower than double exponential in time. It is also interesting to 
examine the degree of nonlinearity in the vortex stretching term during this stage. An 
0{T — t)^^ blowup rate in the maximum vorticity would imply that the nonlinearity in the 
vortex stretching term is quadratic. However, our numerical results show that the vortex 
stretching term, when projected to the unit vorticity vector, is bounded by ||u;||oo log(||cJ||oo), 
where lj is vorticity. It is easy to see that such upper bound on the vortex stretching term 
implies that the maximum vorticity is bounded by double exponential in time. During the 
final stage for 17 < t < 19, we observe that the growth of the maximum vorticity slows down 
considerably and deviates from double exponential growth, indicating that there is stronger 
cancellation taking place in the vortex stretching term. 

We also find that the vortex lines near the region of the maximum vorticity are relatively 
straight and the vorticity vectors seem to be quite regular. This was also observed in . On 
the other hand, the inner region containing the maximum vorticity does not seem to shrink 
to zero at a rate of (T — t)^/^, as predicted by Kerr's computations. Moreover, we find that 
the velocity field, the enstrophy, and enstrophy production rate remain bounded throughout 
the computations. The fact that the velocity field remains bounded is significant. With 
velocity field being bounded, the result of Deng-Hou-Yu |9j can be applied, which implies 
the non-blowup of the Euler equations up to T = 19. The geometric regularity of the vortex 
lines near the inner region seems to play an important role in the dynamic depletion of vortex 
stretching P ITU]. 

We would like to stress the importance of sufficient resolution in determining the nature 
of the nearly singular behavior of the 3D Euler equations. As demonstrated by our numerical 
computations, the 3D Euler equations have different growth rates in the maximum vorticity 
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in different stages. A resolution witliout tlie proper level of refinement would not be able to 
capture the transition from one growth phase to another. In |^ , the inverse of the maximum 
vorticity was shown to approach to zero almost linearly in time up to t = 17. If this trend 
were to continue to hold, it would lead to the blowup of the maximum vorticity in the form 
of 0(T — t)^^. However, with increasing resolutions, we find that the curve corresponding 
to the inverse of maximum vorticity starts to turn away from zero starting at t = 17. We 
also observe that the velocity field becomes saturated around this time. Incidentally, this is 
precisely the time when Kerr's computations began to lose resolution. At t = 17, the thin 
vortex sheets have already started to roll up. After t = 17.5, the vorticity in the rolled up 
region has developed large gradients in all three directions, with z being the most singular 
direction and y being the least singular direction. To resolve the large gradients in all three 
directions, we allocate 3072 grid points along the z direction, 1536 along the x direction, and 
1024 along the y direction. This level of resolution ensures that we have about 8 grid points 
across the most singular region in each direction toward the end of the computations. 

Kerr interpreted the roll-up of the vortex sheet as "two vortex sheets meeting at an angle" 
and argued that the formation of this angle may be responsible for the finite time blowup of 
the Euler equations. Our computations indicate that the roUup region of the vortex sheet is 
still relatively smooth even during the final stage of the computations. Moreover, according 
to the results in ^ ^] , it is the curvature of the vortex lines and the divergence of the 
unit vorticity vector that contribute to the blow-up, not the curvature of the vortex sheet 
itself. Further, we observe that the vortex lines near the region of the maximum vorticity 
are relatively smooth. This geometric regularity leads to strong dynamic depletion of the 
nonlinear vortex stretching. 

There is reason to believe that if the current scenario persists, there is no blowup of 
the 3D Euler equations for these data beyond T = 19. In fact, during the final stage of 
the computations for 17.5 < t < 19, the vortex lines near the region of the maximum 
vorticity remain smooth. Further, as the vortex sheet rolls up, we observe that the location 
of the maximum vorticity moves away from the dividing plane separating the two vortex 
tubes toward the rolled up portion of the vortex sheet, leading to a slower growth rate of 
maximum vorticity. 

The rest of this paper is organized as follows. We describe the set-up of the initial 
condition in Section 2 and describe our numerical method in Section 3. In Section 4, we 
describe our numerical results in detail and perform comparisons with the previous results 
obtained in ^JE|- Some concluding remarks are made in Section 5. 



The 3D incompressible Euler equations in the vorticity stream function formulation are given 



2 The Initial Condition 



as follows: 



UJt + {u ■ V)uj 



Vn ■ uj 



(1) 
(2) 
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with initial condition uj \t=o= cjoj where m is velocity, u is vorticity, and ip is stream function. 
Vorticity is related to velocity by = V x -u. The incompressibility implies that 

V-M = V- d; = V- V^=0. 

We consider periodic boundary conditions with period 27r in all three directions. 

We study the interaction of two perturbed antiparallel vortex tubes using the same initial 
condition as that of Kerr (see Section III of ^1]). There are a few misprints in the analytic 
expression of the initial condition given in jHj. In our computations, we use the corrected 
version of Kerr's initial condition by comparing with Kerr's Fortran subroutine which was 
kindly provided to us by him. A list of corrections to these misprints is given in the Appendix. 

The initial condition is given by a pair of perturbed anti-parallel vortex tubes, which is 
expressed in terms of vorticity. The vorticity that describes the vortex tube above the x-y 
plane is of the form: 

U = Uj{r){uj^, UJy, UJ;,). (3) 

The first step in setting up the initial condition is to define the profile a;(r). If r > 1, we set 
uj{r) = 0. For r < 1, we define 

^(r) =exp[/(r)], (4) 

with /(r) given by 

The radius r is centered around an initial vortex core trajectory (X, Y, Z), and is defined by 

r = \{x,y,z)-{X,Y,Z)\/R, for r ^ 1. (6) 
The initial vortex core trajectory (X, Y, Z) is characterized by 

iX,Y,Z) = ixis),y,zis)), (7) 

where s is a function of y and 

x{s) = Xq + SxCos^ns/Lx), (8) 

z{s) = zq + 6z cos^ns/Lz). (9) 

To complete the definition of Lj{r), we need to define s as a function of y, which is given 
below: 

siy) = 1/2 + LySyi sm{7iy2/Ly), (10) 

and 

y2 = y + Ly5y2 sin {Try/ Ly) . (11) 
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Figure 1: The axial vorticity (uy) contours of the initial value on the symmetry plane. 



The second step is to define the vorticity vector {ujx,ujy,ujz), which is given as follows: 



UJy 



1 + ■KOy.^ COS \ — 

Ly 



1 + TrSy^ COS 



iry 



X 



X 



1 + TTOyi COS I -;— 



1 + TT6yi COS 



sm 



sm 



7Ts{y) 



(12) 
(13) 
(14) 



We choose exactly the same parameters as in [H]. Specifically, we set 5^1 = 0.5, 5y2 = 0.4, 
6x = —1.6, 6z = 0, zq = 1.57 and R = 0.75. The constant xq fixes the center of perturbation 
for the vortex tube along the x direction. In our computations, we set xq = 0. Moreover, we 
choose Lx = Ly = Atc, and Lz = 27i. 

The third step is to rescale the initial profile defined above. According to (see the 
last paragraph on page 1728 of [T¥), we need to rescale the above initial vorticity profile by a 
constant factor so that the maximum vorticity in the y direction is increased to 8. With the 
choice of the above parameters, the maximum vorticity in the y direction before rescaling is 
equal to 0.999766. Thus, the constant rescahng factor is equal to 8.001873. 

The final step in defining the initial condition is to filter the initial vorticity profile. After 
rescaling the initial vorticity, we apply exactly the same Fourier filter as the one used in jl4j . 
i.e. exp [—0.05 (fc^ + ky + kf)~\ , to the initial vorticity to smooth the rough edges. 

We should point out that due to the difference between ours and Kerr's discretization 
strategies in solving the 3D Euler equations, the discrete initial condition generated by Kerr's 
discretization and the one generated by our pseudo-spectral discretization are not exactly the 
same. In Kerr used the Chebyshev polynomials to approximate the solution along the z 
direction. In order to prepare the initial data that can be used for the Chebyshev polynomials, 
he performed some interpolation and used extra filtering. This interpolation and extra 
filtering seem to introduce some asymmetry to Kerr's discrete initial data. According to fl] 
(see the top paragraph of page 1729), "An effect of the initial filter upon the vorticity contours 
at t = is a long tail in Fig. 2(a)". Since we perform pseudo-spectral approximations in 
all three directions, there is no need to perform interpolation or use extra filtering as was 
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Figure 2: The axial vorticity [ujy) contours when t = 6 on the symmetry plane. 

done in [T3]. To demonstrate this slight difference between Kerr's initial data and ours, we 
plot the initial vorticity contours along the symmetry plane in Figure As we can see, the 
initial vorticity contours in Figure Q are essentially symmetric. This is in contrast to the 
apparent asymmetry in Kerr's initial vorticity contours (see Fig. 2(a) of J^). The 3D plot 
of the initial vortex tubes is given in Figure El Again, we can see that the initial vortex 
tube is essentially symmetric. As time increases, the two antiparallel perturbed vortex tubes 
approach each other. By time t = 6, we already observe a significant flattening near the 
center of the tubes, see Figure |21 and Figure El 

3 The Numerical Method 

We use the pseudo-spectral method with a very high order Fourier smoothing to discretize 
the 3D Euler equations. The Fourier smoothing that we use along the Xj direction is of 
the form: p{2kj/Nj) = exp{—a{2kj/Nj)'^) with a = 36 and m = 36, where kj is the wave 
number {\kj\ ^ Nj/2) along the xj direction and Nj is the total number of grid points along 
the Xj direction. Specifically, if % is the discrete Fourier transform of v, then we approximate 
Vxj by taking the discrete inverse Fourier transform of ikjp{2kj/Nj)vk, where k = {ki, k2, k^). 
The time integration is performed using the classical fourth order Runge-Kutta method. 
Adaptive time stepping is used to satisfy the CFL stability condition with CFL number 
equal to 7r/4. We use up to 1536 x 1024 x 3072 space resolution in order to resolve the nearly 
singular behavior of the 3D Euler equations. 

There is a good reason why we choose to use the pseudo-spectral method with the above 
high order Fourier smoothing instead of using the classical 2/3 dealiasing rule. The Fourier 
smoothing we use is designed to keep the majority of the Fourier modes unchanged and 
remove the very high modes to avoid the aliasing errors, see Fig. jHfor the profile of p{x). We 
choose a to be 36 to guarantee that p{k) reaches the level of the round-off errors (O(10~^^)) 
at the highest modes. We choose the order of smoothing, m, to be 36 in order to optimize the 
accuracy of the spectral approximation, while still keeping the aliasing errors under control. 
As we can see from FigurejH the effective modes in our computation are about 12 ~ 15% more 
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Figure 3: The 3D view of the vortex tube for t = and t = 6. The tube is the isosurface 
at 60% of the maximum vorticity. The ribbons on the symmetry plane are the contours at 
other different values. 

than those using the standard 2/3 dealiasing rule. To demonstrate that the pseudo-spectral 
method with the above high order Fourier smoothing is indeed more accurate than the 
pseudo-spectral method with the 2/3 dealiasing rule, we perform resolution study of the two 
approaches. In Figure El we compare the Fourier spectra of the enstrophy obtained by using 
the pseudo-spectral method with the 2/3 dealiasing rule with those obtained by the pseudo- 
spectral method with the high order smoothing. For a fixed resolution 768 x 512 x 1536, 
we can see that the Fourier spectra obtained by the pseudo-spectral method with the high 
order smoothing are more accurate than those obtained by the spectral method using the 
2/3 dealiasing rule. This can be seen by comparing the results with the corresponding 
computations using a higher resolution 1024 x 768 x 2048. Moreover, the pseudo-spectral 
method using the high order Fourier smoothing does not give the spurious oscillations in the 
Fourier spectra which are present in the computations using the 2/3 dealiasing rule near the 
2/3 cut-off point. 

We have used a sequence of resolutions in our computations in order to perform a resolu- 
tion study. The resolutions presented in the paper are 768 x 512 x 1536, 1024 x 768 x 2048, 
and 1536 x 1024 x 3072 respectively. Except for the computation on the largest resolution 
1536 X 1024 X 3072, all computations are carried out from t = to t = 19. The computation 
on the final resolution 1536 x 1024 x 3072 is started from t = 10 with the initial condition 
given by the computation with the resolution 1024 x 768 x 2048. From our resolution study 
for the velocity and vorticity in both physical and spectral spaces, we find that the solution 
at if: = 10 is fully resolved even on the resolution 768 x 512 x 1536. Thus it is justified to 
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Figure 4: The profile of tlie Fourier smootliing: exp(— 36(x)'^*^). Tlie vertical line corre- 
sponds to the cut-off mode using the 2/3 dealiasing rule. We can sec that using this Fourier 
smoothing we keep about 12 ~ 15% more modes than those using the 2/3 deahasing rule. 

start the computation with the largest resolution at i = 10 using the computation obtained 
by the resolution 1024 x 768 x 2048 as the initial condition. 

We also exploit the symmetry properties of the solution in our computations, and perform 
our computations on only a quarter of the whole domain. We use the well-known parallel 
FFT package, called FFTW 3.1 to perform the sine and cosine transformations. The whole 
program is coded in C and we use MPI as the parallel interface. Since the solution appears 
to be most singular in the z direction, we allocate twice as many grid points along the z 
direction than along the x direction. The solution is least singular in the y direction. We 
allocate the smallest resolution in the y direction to reduce the computation cost. In our 
computations, two typical ratios in the resolution along the x, y and z directions are 3:2:6 
and 4 : 3 : 8. 

In the early stage of the computations for < t < 8, the solution is still very regular. 
The solution in this stage can be fully resolved using the resolution 768 x 512 x 1536. In 
the intermediate stage for 8 < t < 16, the solution has a fast growth in the maximum 
vorticity. After t = 16, the two vortex tubes become severely flattened, and evolve into two 
thin vortex sheets, which roll up subsequently. The maximum vorticity also experiences a 
transition from the double exponential growth to a slower growth rate due to the dynamic 
depletion of vortex stretching. In order to resolve the later stage of the solution, higher 
resolutions become necessary. 

Our computations were carried out on the PC chistcr LSSC-II in the Institute of Compu- 
tational Mathematics and Scientific/Engineering Computing of Chinese Academy of Sciences 
and the Shentcng 6800 cluster in the Super Computing Center of Chinese Academy of Sci- 
ences. The maximal memory consumption in our computations is about 120 GBytes. 
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Figure 5: The comparison of the enstrophy spectra obtained using the high order Fourier 
smoothing method with those using the 2/3 deahasing rule. The dashed hues and dashed- 
dotted hues are the enstrophy spectra with the resolution 768 x 512 x 1536 using the 2/3 
dealiasing rule and the Fourier smoothing, respectively. The solid lines are the enstrophy 
spectra with resolution 1024 x 768 x 2048. The times for the spectra lines are at t = 
15, 16, 17, 18, 19 respectively. 

4 Numerical Results 
4.1 Review of Kerr's results 

In pi], Kerr presented numerical evidence which suggested a finite time singularity of the 
3D Euler equations for two perturbed antiparallel vortex tubes. He used a pseudo-spectral 
discretization in the x and y directions, and a Chebyshev method in the z direction with 
resolution of order 512 x 256 x 192. His computations showed that the growth of the peak 
vorticity, the peak axial strain, and the enstrophy production obey (T — t)~^ with T = 
18.9. As t approaches the alleged blow-up time T, the region bounded by the contour 
of 0.6||a;||oo, also known as the inner region, looks like two vortex sheets with thickness 
~ (T — t) meeting at an angle ^3]. This region has the length scale (T — t)^/^ in the 
vorticity direction. The maximum vorticity resides in the small tube-like region, with scaling 
(T — t)^/^ X {T — t) X (T — t), which is the intersection of the two sheets. Inside the inner 
region, the vortex lines are "relatively straight" J3] . Kerr stated in his paper ^1] (see page 
1727) that his numerical results shown after t = 17 and up to t = 18 were "not part of the 
primary evidence for a singularity" due to the lack of sufficient numerical resolution and the 
presence of noise in the numerical solutions. 

In his recent paper jTT] (see also fSlEl); Kerr applied a high wave number filter to the 
data obtained in his original computations to "remove the noise that masked the structures 
in earlier graphics" presented in ^3]. With this filtered solution, he presented some scaling 
analysis of the numerical solutions up to t = 17.5. Two new properties were presented in 
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this recent paper ^7j. First, the velocity field was shown to blow up like 0{T — t)"^/^ with 
T being revised to T = 18.7. The maximum velocity is located on the boundary of the 
inner region with distance (T — t)^/^ away from the position where the maximum vorticity 
is achieved. Secondly, he showed that the blowup is characterized by two anisotropic length 
scales, p ~ (T — t) and i? ~ (T — t)^/^. Further, Kerr argued that the curvature, n, of the 
vortex lines and V ■ ^ (here ^ = u;/|a;|) in the inner region are likely bounded by (T — t)^^/^ 



4.2 Recent theoretical results on the 3D Euler equations 

There has been some interesting development in the theoretical understanding of the 3D 
incompressible Euler equations. It has been shown by several authors that the local geometric 
regularity of vortex lines can play an important role in depleting nonlinear vortex stretching. 
In particular, Constantin, Fefferman and Majda [HJ proved that if (i) there is up to time 
T an 0(1) region Q in which the vorticity vector is smoothly directed, i.e., the maximum 
norm of in this 0(1) region is integrable in time, (ii) the maximum norm of velocity 
is uniformly bounded in time, plus a technical condition on the distribution of the vorticity 
within this region, then no blow-up can occur in this region up to time T. While this result 
is very interesting, it does not apply to Kerr's computations since the two assumptions are 
violated by Kerr's computations. 

Motivated by the result of [Q, Deng, Hou, and Yu [9|| obtained sharper non-blowup 
conditions which use only very localized information of the vortex lines. Assume that at 
each time t there exists some vortex line segment Lf on which the local maximum vorticity 
is comparable to the global maximum vorticity. Further, we denote L{t) as the arclength 
of Lf. Roughly speaking, if (i) the velocity field along Lt is bounded by Cu(T — t)~" for 
some a < 1, (ii) Cl{T — tY < L{t) < Cq/ max^^d/tl, |V ■ .^|), for some P < 1 — a, then the 
solution of the 3D Euler equations remains regular up to T. In Kerr's computations, the first 
condition is satisfied with a = 1/2. Moreover, based on the bound on k, and V-^ in the inner 
region ^Tj, we can choose a vortex line segment of length (T — tY^"^ in the inner region so that 
the upper bound in the second condition is satisfied. However, the lower bound is violated 
since (3 < 1/2 = 1 — a. In a subsequent paper the non-blowup conditions have been 
improved to include the critical case, = 1 — a by requiring the scahng constants, Cu, Cq, 
and Oj;, in conditions (i)-(ii), to satisfy an algebraic inequality. This algebraic inequality can 
be checked numerically if we obtain a good estimate of these scaling constants. For example, 
if Co = 0.1, which seems reasonable since the vortex lines are relatively straight in the inner 
region, the result of ITU] implies no blowup up to T if 2Cu < 0.430^. However, such scahng 
constants are not available in p^. One of our original motivations was to repeat Kerr's 
computations with higher resolution to obtain a good estimate of these scaling constants. 



11 



4.3 Maximum vorticity growth 

We first present the result on the growth of the maximum vorticity in time, see Figure IHl 
The maximum vorticity increases rapidly from the initial value of 0.669 to 23.46 at the final 
time t = 19, a factor of 35 increase from its initial value. Kerr's computations predicted a 
finite time singularity at T = 18.7. Our computations show no sign of finite time blowup 
of the 3D Euler equations up to T = 19, beyond the singularity time predicted by Kerr. 
In Figure El we plot the maximum vorticity in time using three different resolutions, i.e. 
768 X 512 X 1536, 1024 x 768 x 2048, and 1536 x 1024 x 3072 respectively. As we can see, the 
agreement between the two successive resolutions is very good with only mild disagreement 
toward the end of the computations. This indicates that a very high space resolution is 
indeed needed to capture the rapid growth of maximum vorticity at the later stage of the 
computations. 

We observe that the growth of the maximum vorticity has three distinguished phases. 
The first stage is for < t < 12. In this early stage, the maximum vorticity grows only 
exponentially in time. This is consistent with Kerr's results. The second stage is for 12 < 
t < 17. During this intermediate stage, the two vortex tubes experience tremendous core 
deformation and become severely flattened. Each vortex tube effectively turns into a vortex 
sheet with rapidly decreasing thickness. We observe that the growth of maximum vorticity 
is slightly slower than double exponential in time during the second stage, see Figure [7| 
This growth behavior can be also confirmed by examining the degree of nonlinearity in the 
vortex stretching term. If the maximum vorticity indeed blew up like 0{T — t)~^, as alleged 
in ^3], the vortex stretching term at the position of the maximum vorticity should have 
been quadratic as a function of maximum vorticity. However, as Figure |H1 shows, the vortex 
stretching term, when projected to the unit vorticity vector, grows much slower than the 
quadratic nonlinearity. In fact, it is even slower than C||u;||oo log(||a;||oo), i.e. 

lie- Vm-cD|U < C7p|Ulog(p||oo), 15<t<19. (15) 
Using the equation that governs the magnitude of vorticity [7], 

d 

— \u\ + {u-\/)\u\=^-Vu-Cu, (16) 

one can easily show that inequality (fTH|l implies that the maximum vorticity is bounded by 
double exponential in time. 

During the final stage for 17 < t < 19, we observe that the growth of the maximum vor- 
ticity slows down and deviates from double exponential growth, see Figure [7| This indicates 
that there is stronger cancellation taking place in the vortex stretching term. 

We remark that for vorticity that grows as rapidly as double exponential in time, one 
may be tempted to fit the maximum vorticity growth by (T — t)^^ for some T. Indeed, if we 
choose T = 18.7 as suggested by Kerr in ^Jj, we find a reasonably good fit for the maximum 
vorticity as a function of c/ (T — t) for the period 15 < t < 17. We plot the scaling constant 
c in Figure ini As we can see, c is close to a constant for 15 < t < 17. To conclude that the 
3D Euler equations indeed develop a finite time singularity, one must demonstrate that such 
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Figure 6: The maximum vorticity ||ct;||oo in time using different resolutions. 




Figure 7: The plot of log log ||ct;||oo vs time, resolution 1536 x 1024 x 3072. 
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Figure 8: Study of the vortex stretching term in time, resolution 1536 x 1024 x 3072. We 
take Ci = 1/8.128, C2 = 1/23.24 to match the same starting value for all three plots. 

scaling persists as t approaches to T. As we can see from Figure IHl the scaling constant c 
decreases rapidly to zero as t approaches to the alleged singularity time T. Therefore, the 
fitting of ||u;||oo ~ 0{T — t)~^ is not correct asymptotically. 

A similar test can be performed for the inverse of the maximum vorticity. In Figure 
inn we plot the inverse of the maximum vorticity using different resolutions. As we can see 
from this picture, the inverse of the maximum vorticity approaches to zero almost linearly 
in time for 8 < t < 17. This was one of the strong evidences presented in that suggests 
a finite time blowup of the 3D Euler equations. If this trend were to continue to hold up 
to T, it would have led to the blowup of the maximum vorticity in the form of 0{T — t)^^. 
However, as we increase our resolutions, we find that the curve corresponding to the inverse 
of the maximum vorticity starts to turn away from zero around t = 17. This is precisely 
the time when Kerr's computations began to lose resolution. By t = 17.5, the gradients of 
the solution become very large in all three directions. In order to resolve the nearly singular 
solution structure, we use 1536 x 1024 x 3072 grid points from t = 10 to 19. This level of 
resolution gives about 8 grid points across the most singular region in each direction toward 
the end of the computations. 

4.4 Velocity profile 

One of the important findings of our computations is that the velocity field is actually 
bounded by 1/2 up to T = 19. This is in contrast to Kerr's computations in which the 
maximum velocity was shown to blow up like 0{T — We plot the maximum velocity 
as a function of time using different resolutions in Figure ^2 The computation with the 
resolution 768 x 512 x 1536 shows some mild discrepancy toward the end of the computation. 
On the other hand, the computation obtained by resolution 1024 x 768 x 2048 and the one 
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Figure 9: Scaling constant in time for the fitting ||a;||oo ~ c/{T — t), T — 18.7. 
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Figure 10: Tlie inverse of maximum vorticity ||(^||oo in time using different resolutions. 
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Figure 11: Maximum velocity \\u\\^ in time using different resolutions. 

obtained by resolution 1536 x 1024 x 3072 are almost indistinguishable. As we can see, the 
maximum velocity grows slowly in time and is relatively small in magnitude. There is a 
relatively fast growth of maximum velocity between t = 14 and t = 17. But this growth 
becomes saturated hj t = 17.5. After t = 18.4, the velocity experiences a mild growth, but 
it is still bounded by 0.46 at the final time T = 19. We also plot the contours of \u\ near 
the region of maximum vorticity at t = 18 and 19 in Figure IT^ As we can see, the velocity 
seems to be well resolved. 

The fact that the velocity field is bounded is significant. By re-examining the non- 
blowup conditions of the theory of Deng-Hou-Yu [HI, we find that the first condition is 
now satisfied with a = since the velocity field is bounded. According to ^7j, we have 
max^jd^l, |V ■ ^1) < Co(T — 1)^^^"^. In fact, our computations indicate that the curvature 
and the divergence of the unit vorticity vector are actually bounded. With a = 0, we can 
now choose the vortex line segment of length, L{t) = (T — t)^ with /3=1/2<1 — a, so 
that the second condition is now satisfied. Thus, the theory of Deng-Hou-Yu applies, which 
imphes non-blowup of the 3D Euler equations up to T. 

4.5 Local vorticity structure 

In this subsection, we would like to examine the local vorticity structure near the region 
of the maximum vorticity. To illustrate the development in the symmetry plane, we show 
a series of vorticity contours near the region of the maximum vorticity at late times in a 
manner similar to the results presented in jHj. In Kerr's computations, he observed that 
the head and tail in the symmetry plane develop a corner separating the head and tail. We 
adopt Kerr's definition of the "head" to be the region extending above the vorticity peak Up 
just behind the leading edge of the vortex sheet. The "tail" is the vortex sheet extending 
behind the peak vorticity. One interesting question is to determine whether one direction 
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Figure 12: The contour of 1^1 in the region of maximum vorticity on the symmetry plane. 
The pictures are plotted at t = 18 and t = 19 respectively using resolution 1536 x 1024 x 3072. 



becomes progressively more flattened or stretched as the flow evolves and whether the rates 
of collapse are the same in different directions. Our computational results are in qualitative 
agreement with Kerr's in the early and intermediate stages. In particular, we observe that 
as the flow evolves the region of peak vorticity concentrates into the region where the vortex 
sheets of the head and tail meet. To compare with Kerr's figures, we scale the vorticity 
contours in the x — z plane by a factor of 5 in the z direction. The results at t = 15 and 
t = 17 are plotted in Figure El We can see that the location of maximum axial vorticity 
moves toward the corner where the vortex sheets of the head and tail meet as time increases, 
see also Figure IT^ This is in qualitative agreement with Kerr's results. 

In order to see better the dynamic development of the local vortex structure, we plot a 
sequence of vorticity contours on the symmetry plane at t = 16, 17, 18, and 19 respectively 
in Figure HH The pictures are plotted using the original length scales, without the scaling 
by a factor of 5 in the z direction as in Figure El From these results, we can see that 
the vortex sheet is compressed in the z direction. It is clear that a thin layer (or a vortex 
sheet) is formed dynamically. The head of the vortex sheet is a bit thicker than the tail at 
the beginning. The head of the vortex sheet begins to roll up around t = 16. By the time 
t = 19, the head of the vortex sheet has traveled backward for quite a distance, and the 
vortex sheet has been compressed quite strongly along the z direction. 

We also plot the isosurface of vorticity near the region of the maximum vorticity in 
Figures El and El to illustrate the dynamic roll-up of the vortex sheet near the region of 
the maximum vorticity. Figure El gives the local vorticity structure at t = 17. If we scale 
the local roll-up region on the left hand side next to the box by a factor of 4 along the z 
direction, as was done in ^\ , we would obtain a local roll- up structure which is qualitatively 
similar to Figure 1 in In Figure El we show the local vorticity structure for t = 18 and 
19. In both figures, the isosurface is set at 0.5 x ||cJ||oo. Here we make a few observations. 
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Figure 13: The contour of axial vorticity around the maximum vorticity on the symmetry 
plane at t = 15, 17. The figure is scaled in z direction by a factor of 5 to compare with the 
Figure 4 in [14j. A contour at value very close to the maximum value is plotted to show the 
location of the maximum vorticity. 




Figure 14: The contour of axial vorticity around the maximum vorticity on the symmetry 
plane at t = 17.5, 18, 18.5, 19. 
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Figure 15: The local 3D vortex structure and vortex lines around the maximum vorticity at 
t = 17. The size of the box on the left is 0.075'^ to demonstrate the scale of the picture. 

First, the vortex lines near the region of maximum vorticity are relatively straight and the 
vorticity vectors seem to be quite regular. This was also observed in ^3]. On the other 
hand, the inner region containing the maximum vorticity does not seem to shrink to zero at 
a rate of (T — t)^/^, as predicted in The length and the width of the vortex sheet are 
still 0(1), although the thickness of the vortex sheet becomes quite small. 

We also plot the energy spectrum in Figure[T7|at t = 16, 17, 18, 19. A finite time blow-up 
of enstrophy would imply that the energy spectrum decays no faster than \k\~^. In [H], 
Kerr observed that the energy spectrum decays exactly like \k\~^, suggesting a finite time 
blow-up of the enstrophy (recall that enstrophy is defined as the square of the norm of 
vorticity, i.e. ||a;|||). Our computations show that the energy spectrum approaches \k\~^ for 
\k\ < 100 as time increases to t = 19. This is in qualitative agreement with Kerr's results. 
Note that there are only less than 100 modes available along the \kx\ or \ky\ direction in 
Kerr's computations, see Figure 18 (a)-(b) of ^j. On the other hand, our computations 
show that the high frequency Fourier spectrum for 100 < < 1300 decays much faster 
than \k\~^, as one can see from Figures IT7I and 1^ This indicates that there is no blow-up 
in enstrophy. This is also supported by the enstrophy spectrum, given in Figure ^1 and the 
plot of enstrophy as a function of time in Figures ^1 In Figure 1201 we plot the enstrophy 
production rate, which is defined as the time derivative of enstrophy. Although it grows 
relatively fast, it actually grows slower than double exponential in time (see the picture on 
the right in Figure EOI)- In the double logarithm plot of the enstrophy production rate, we 
multiply the enstrophy production rate by a constant factor 8 to make the second logarithm 
well-defined. It is interesting to note that the double logarithm of enstrophy production rate 
in Figure QUI is qualitatively similar to the double logarithm of ||ci3||oo in Figure [3 
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Figure 16: The local 3D vortex structures and vortex lines around the maximum vorticity 
at t = 18 (on the left) and t = 19 (on the right). 




Figure 17: The energy spectra for velocity at different times in log-log scale. The energy 
spectrum is calculated using E{k) = J2\k\£{k-i/2 k+1/2] l^fcl^- "^^^ time for the spectral lines 
from bottom to top are t = 15, 16, 17, 18, 19. The dashed line corresponds to k~^. 
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Figure 18: The enstrophy spectra for vorticity at different times in log-log scale. The en- 
strophy spectrum is calculated using Ens{k) — 'Yli\k\^{k-i/2k+i/2]\^k\^- '^^^ dashed line 
corresponds to k~^. 
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Figure 19: The enstrophy as a function of time, resolution 1536 x 1024 x 3072. 
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Figure 20: The enstrophy production rate (left) and the double logarithm of 8 times enstro- 
phy production rate (right) as a function of time, resolution 1536 x 1024 x 3072. 

4.6 Resolution Study 

In this subsection, we perform a resolution study to make sure that the nearly singular 
behavior of the 3D Euler equations is resolved by our computational grid. In Figure El we 
have performed resolution study for the maximum vorticity using three different resolutions 
and found very good agreement for the time interval [0, 18]. There is only a mild disagreement 
toward the end of the computations from t = 18 to 19. In our computations with different 
resolutions, we find that the maximum vorticity always grows slower than double exponential 
in time. Similar resolution study has been performed for the inverse of maximum vorticity 
in Figure ^1 and for the maximum velocity in Figure ^2 We observe excellent agreement 
between the solutions obtained by the two largest resolutions. We have also performed 
similar resolution study in the Fourier space by examining the convergence of the energy 
and enstrophy spectra using different resolutions. We observe that the Fourier spectrum 
corresponding to the effective modes of one resolution is in excellent agreement with that 
corresponding a higher resolution computation, see Figures 1221 and |221 

To see how many grid points we have across the most singular region, we plot the under- 
lying mesh for the vorticity contours in the x — z plane in Figure |^ One can see from this 
picture that we have about 16 grid points in the z direction at t = 18 and 8 grid points at 
t = 19. It is also interesting to note that at t = 18.5, the location of the maximum vorticity 
has moved away from the bottom of the vortex sheet structure. If the current trend contin- 
ues, it is likely that the location of the maximum vorticity will continue to move away from 
the bottom of the vortex sheet. One of the possible blow-up scenarios is that the interaction 
of the two perturbed antiparallel vortex tubes would induce a strong compression between 
the two vortex tubes, leading to a finite time collapse of the two vortex tubes. The fact that 
the location of the maximum vorticity moves away from the dividing plane of the two vortex 
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Figure 21: This picture is to illustrate the mesh around the maximum vorticity. The times 
for this plot are 18, 19. At t = 19, we still have about 8 points along z direction to resolve 
the nearly singular layer. 

tubes seems to destroy the desired mechanism to produce a blow-up. 

5 Concluding Remarks 

We investigate the interaction of two perturbed vortex tubes for the 3D Euler equations 
using Kerr's initial data. Our numerical computations demonstrate a very subtle dynamic 
depletion of vortex stretching. The maximum vorticity is shown to grow no faster than 
double exponential in time up to T = 19, beyond the singularity time predicted by Kerr in 
[Tij . The local geometric regularity of vortex lines seems to be responsible for this dynamic 
depletion of vortex stretching. Sufficient numerical resolution is essential in capturing the 
double exponential growth in vorticity and the dynamic depletion of vortex stretching. The 
velocity field and the enstrophy are shown to be bounded throughout the computations. We 
provide evidence that the vortex stretching term is only weakly nonlinear and is bounded 
by ||u;||oo log(||w||oo)- Such an upper bound on the vortex stretching term implies that the 
maximum vorticity is bounded by the double exponential in time. Our computational results 
also satisfy the non-blowup conditions of Deng-Hou-Yu, which provides a theoretical support 
for our computational results. 

The current computations, even with this level of resolution, can not rule out the possi- 
bility of the blow-up of the 3D Euler equations for large times for Kerr's initial data. The 
theoretical results of jHl IHl UHj and the computations presented here suggest that a finite time 
singularity, if it exists, would have rather complicated geometric structures. There are other 
types of potential Euler singularities that are not considered in this paper. Among them. 
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Figure 22: Convergence study for enstrophy spectra using different resolutions. The dashed 
hues and the sohd hues are the enstrophy spectra on resolution 1536 x 1024 x 3072 and 1024 x 
768 x 2048, respectively. The times for the hnes from bottom to top are t — 16, 17, 18, 19. 




Figure 23: Convergence study for energy spectra using different resolutions. The dashed 
hnes and the sohd lines are the energy spectra on resolution 1536 x 1024 x 3072 and 1024 x 
768 X 2048, respectively. The times for the hnes from bottom to top are t — 16, 17, 18, 19. 
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the Kida-Pelz initial condition [3 1^ is worth further investigation. The extra symmetry 
constraints in this type of initial data are believed to be important in producing a finite time 
singularity for the 3D Euler equations. Indeed, the computations by Boratav and Pelz |3] 
and Pelz [21] indicate a more singular self-similar type of blow-up. Pelz's computations also 
fall in the critical case of the non-blowup theory of Deng-Hou-Yu We are currently 

investigating this problem numerically using even higher resolutions. We will report the 
results elsewhere. 

Appendix. Corrections to Some Misprints in [114] 

In this appendix, we explain the corrections that we make regarding the misprints in the 
description of the initial condition in jT3] . There are two constraints on the initial vorticity. 
The first one is that it must be divergence free. The second one is that it must satisfy the 
periodic boundary condition. It is obvious that we have 

V ■ {ux, Uy, u^) = 0. (17) 

Thus the divergence free constraint on the initial vorticity implies that uj{r) must satisfy 

Vuj{r) ■ {u^, Uy, u;^) = 0. 

The analytic expression of the initial vorticity profile in ^3] does not satisfy the above 
constraint due to a few typos in various formula. We correct these typos by comparing the 
analytic formula with the formula that were actually used by Kerr in his Fortran subroutine 
that generates the initial data. Below we would list these typos and point out the corrections. 

1. In equation ((Tj), the original expression in [T3] was written as [xq + x{s),y, Zq + z{s)]. 
We remove xq and zq from ((Tj) since the definition of x{s) and z{s) has already taken 
Xq and Zq into account. 

2. In (jHJ, the original expression in |14j was given as Xq + (5^cos(s). This would violate 
the divergence free condition. We correct this by replacing cos(s) with cos (tcs/Lx). 

3. In (jH)), the original expression was zq + 6zCos{s). This again violates the divergence 
free condition. We correct it by replacing cos(s) with cos{tts/Lz). 

4. In ()12|). the last factor in the original expression was sin(7rs(?/)). We correct it by 
replacing sin(7rs(?/)) with sin {ns^y)/ L^). 

5. In ()14|). the last factor in the original expression was sin(7rs(?/)). We correct it by 
replacing sm{ns{y)) with sin {ns{y)/ L^). 

With the above corrections, it can be verified that both the divergence free condition and 
the periodic boundary conditions are satisfied. 

Based on Kerr's subroutine that generates the initial data, we also make one minor 
modification in the definition of f{r). The original equation in ^3] for (0) was given by 

2 

/(r) = + (1 + + r^) . (18) 
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After studying Kerr's code, we found that the function /(r) in his code was actually defined 
using (0) instead of the above formula. The difference lies in the first factor in the second 
term of the above equation. What was used in Kerr's code was instead in the above 
equation. This minor modification has little effect on the behavior of the solution from 
our computational experience. We make this minor modification to /(r) in order to match 
exactly the initial condition that was actually used in Kerr's computations. 
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